Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae.

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Cyclization mechanism of monoterpenes catalyzed by monoterpene synthases in dipterocarpaceae.

Monoterpenoids are typically present in the secretory tissues of higher plants, and their biosynthesis is catalyzed by the action of monoterpene synthases (MTSs). However, the knowledge about these enzymes is restricted in a few plant species. MTSs are responsible for the complex cyclization of monoterpene precursors, resulting in the production of diverse monoterpene products. These enzymatic reactions are considered exceptionally complex in nature. Therefore, it is crucial to understand the catalytic mechanism of MTSs to elucidate their ability to produce diverse or specific monoterpenoid products. In our study, we analyzed thirteen genomes of Dipterocarpaceae and identified 38 MTSs that generate a variety of monoterpene products. By focusing on four MTSs with different product spectra and analyzing the formation mechanism of acyclic, monocyclic and bicyclic products in MTSs, we observed that even a single amino acid mutation can change the specificity and diversity of MTS products, which is due to the synergistic effect between the shape of the active cavity and the stabilization of carbon-positive intermediates that the mutation changing. Notably, residues N340, I448, and phosphoric acid groups were found to be significant contributors to the stabilization of intermediate terpinyl and pinene cations. Alterations in these residues, either directly or indirectly, can impact the synthesis of single monoterpenes or their mixtures. By revealing the role of key residues in the catalytic process and establishing the interaction model between specific residues and complex monoterpenes in MTSs, it will be possible to reasonably design and engineer different catalytic activities into existing MTSs, laying a foundation for the artificial design and industrial application of MTSs.

Chromosome-level genomes of three key Allium crops and their trait evolution.

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Thirteen Dipterocarpoideae genomes provide insights into their evolution and borneol biosynthesis.

Dipterocarpoideae, the largest subfamily of the Dipterocarpaceae, is a dominant component of Southeast Asian rainforests and is widely used as a source of wood, damar resin, medicine, and essential oil. However, many Dipterocarpoideae species are currently on the IUCN Red List owing to severe degradation of their habitats under global climate change and human disturbance. Genetic information regarding these taxa has only recently been reported with the sequencing of four Dipterocarp genomes, providing clues to the function and evolution of these species. Here, we report on 13 high-quality Dipterocarpoideae genome assemblies, ranging in size from 302.6 to 494.8 Mb and representing the five most species-rich genera in Dipterocarpoideae. Molecular dating analyses support the Western Gondwanaland origin of Dipterocarpaceae. Based on evolutionary analysis, we propose a three-step chromosome evolution scenario to describe the karyotypic evolution from an ancestor with six chromosomes to present-day species with 11 and 7 chromosomes. We discovered an expansion of genes encoding cellulose synthase (CesA), which is essential for cellulose biosynthesis and secondary cell-wall formation. We functionally identified five bornyl diphosphate synthase (BPPS) genes, which specifically catalyze the biosynthesis of borneol, a natural medicinal compound extracted from damar resin and oils, thus providing a basis for large-scale production of natural borneol in vitro.

Comparison of ONT and CCS sequencing technologies on the polyploid genome of a medicinal plant showed that high error rate of ONT reads are not suitable for self-correction.

BACKGROUND: Many medicinal plants are known for their complex genomes with high ploidy, heterozygosity, and repetitive content which pose severe challenges for genome sequencing of those species. Long reads from Oxford nanopore sequencing technology (ONT) or Pacific Biosciences Single Molecule, Real-Time (SMRT) sequencing offer great advantages in de novo genome assembly, especially for complex genomes with high heterozygosity and repetitive content. Currently, multiple allotetraploid species have sequenced their genomes by long-read sequencing. However, we found that a considerable proportion of these genomes (7.9% on average, maximum 23.7%) could not be covered by NGS (Next Generation Sequencing) reads (uncovered region by NGS reads, UCR) suggesting the questionable and low-quality of those area or genomic areas that can't be sequenced by NGS due to sequencing bias. The underlying causes of those UCR in the genome assembly and solutions to this problem have never been studied. METHODS: In the study, we sequenced the tetraploid genome of Veratrum dahuricum (Turcz.) O. Loes (VDL), a Chinese medicinal plant, with ONT platform and assembled the genome with three strategies in parallel. We compared the qualities, coverage, and heterozygosity of the three ONT assemblies with another released assembly of the same individual using reads from PacBio circular consensus sequencing (CCS) technology, to explore the cause of the UCR. RESULTS: By mapping the NGS reads against the three ONT assemblies and the CCS assembly, we found that the coverage of those ONT assemblies by NGS reads ranged from 49.15 to 76.31%, much smaller than that of the CCS assembly (99.53%). And alignment between ONT assemblies and CCS assembly showed that most UCR can be aligned with CCS assembly. So, we conclude that the UCRs in ONT assembly are low-quality sequences with a high error rate that can't be aligned with short reads, rather than genomic regions that can't be sequenced by NGS. Further comparison among the intermediate versions of ONT assemblies showed that the most probable origin of those errors is a combination of artificial errors introduced by "self-correction" and initial sequencing error in long reads. We also found that polishing the ONT assembly with CCS reads can correct those errors efficiently. CONCLUSIONS: Through analyzing genome features and reads alignment, we have found the causes for the high proportion of UCR in ONT assembly of VDL are sequencing errors and additional errors introduced by self-correction. The high error rates of ONT-raw reads make them not suitable for self-correction prior to allotetraploid genome assembly, as the self-correction will introduce artificial errors to > 5% of the UCR sequences. We suggest high-precision CCS reads be used to polish the assembly to correct those errors effectively for polyploid genomes.

Microbial spore genetic marker technology, a potential technology for traditional Chinese medicine traceability system.

Traditional Chinese medicine (TCM) has a long history, rich clinical experience, and unique advantages in the prevention and treatment of diseases. The quality and safety of Chinese medicinal materials (CMMs) directly affect the clinical efficacy and development of the TCM industry. However, confused provenance, counterfeiting and adulteration of CMMs hinder the acceptance of its therapeutic benefits in modern society. Therefore, the establishment and improvement of a TCM traceability system would be conducive to the transparency of the CMMs production, distribution, and circulation, thereby improving drug safety and promoting industry development. This review discusses the challenges faced in the development of TCM traceability system, the technologies currently available for tracing CMMs, and the potential application of Barcoded Microbial Spores (BMS) to improve CMMs origin traceability and TCM traceability systems.

Chromosome-level genome assembly of Welwitschia mirabilis, a unique Namib Desert species.

Welwitschia mirabilis, which is endemic to the Namib Desert, is the only living species within the family Welwitschiaceae. This species has an extremely long lifespan of up to 2,000 years and bears a single pair of opposite leaves that persist whilst alive. However, the underlying genetic mechanisms and evolution of the species remain poorly elucidated. Here, we report on a chromosome-level genome assembly for W. mirabilis, with a 6.30-Gb genome sequence and contig N50 of 27.50 Mb. In total, 39,019 protein-coding genes were predicted from the genome. Two brassinosteroid-related genes (BRI1 and CYCD3), key regulators of cell division and elongation, were strongly selected in W. mirabilis and may contribute to their long ever-growing leaves. Furthermore, 29 gene families in the mitogen-activated protein kinase signalling pathway showed significant expansion, which may contribute to the desert adaptations of the plant. Three positively selected genes (EHMT1, EIF4E, SOD2) may be involved in the mechanisms leading to long lifespan. Based on molecular clock dating and fossil calibrations, the divergence time of W. mirabilis and Gnetum montanum was estimated at ~123.5 million years ago. Reconstruction of population dynamics from genome data coincided well with the aridification of the Namib Desert. The genome sequence detailed in the current study provides insight into the evolution of W. mirabilis and should be an important resource for further study on gnetophyte and gymnosperm evolution.

Comparative and phylogenetic analyses of eleven complete chloroplast genomes of Dipterocarpoideae.

BACKGROUND: In South-east Asia, Dipterocarpoideae is predominant in most mature forest communities, comprising around 20% of all trees. As large quantity and high quality wood are produced in many species, Dipterocarpoideae plants are the most important and valuable source in the timber market. The d-borneol is one of the essential oil components from Dipterocarpoideae (for example, Dryobalanops aromatica or Dipterocarpus turbinatus) and it is also an important traditional Chinese medicine (TCM) formulation known as "Bingpian" in Chinese, with antibacterial, analgesic and anti-inflammatory effects and can enhance anticancer efficiency. METHODS: In this study, we analyzed 20 chloroplast (cp) genomes characteristics of Dipterocarpoideae, including eleven newly reported genomes and nine cp genomes previously published elsewhere, then we explored the chloroplast genomic features, inverted repeats contraction and expansion, codon usage, amino acid frequency, the repeat sequences and selective pressure analyses. At last, we constructed phylogenetic relationships of Dipterocarpoideae and found the potential barcoding loci. RESULTS: The cp genome of this subfamily has a typical quadripartite structure and maintains a high degree of consistency among species. There were slightly more tandem repeats in cp genomes of Dipterocarpus and Vatica, and the psbH gene was subjected to positive selection in the common ancestor of all the 20 species of Dipterocarpoideae compared with three outgroups. Phylogenetic tree showed that genus Shorea was not a monophyletic group, some Shorea species and genus Parashorea are placed in one clade. In addition, the rpoC2 gene can be used as a potential marker to achieve accurate and rapid species identification in subfamily Dipterocarpoideae. CONCLUSIONS: Dipterocarpoideae had similar cp genomic features and psbM, rbcL, psbH may function in the growth of Dipterocarpoideae. Phylogenetic analysis suggested new taxon treatment is needed for this subfamily indentification. In addition, rpoC2 is potential to be a barcoding gene to TCM distinguish.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau.

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

Development of SSR markers for a Tibetan medicinal plant, Lancea tibetica (Phrymaceae), based on RAD sequencing.

PREMISE OF THE STUDY: Lancea tibetica (Phrymaceae), a Tibetan medicinal plant, is endemic to the Qinghai-Tibet Plateau. The over-exploitation of wild L. tibetica has led to the destruction of many populations. To enhance protection and management, biological research, especially population genetic studies, should be carried out on L. tibetica. Simple sequence repeat (SSR) markers of L. tibetica were developed to analyze population diversity. METHODS AND RESULTS: Four thousand four hundred and forty-one SSR loci were identified for L. tibetica based on restriction-site associated DNA (RAD) sequencing on the Illumina HiSeq platform. One hundred SSR loci were arbitrarily selected for primer design, and 38 of them were successfully amplified. These markers were tested on 56 individuals from three populations of L. tibetica, and 10 markers displayed polymorphisms. The total number of alleles per locus ranged from three to eight, and observed and expected heterozygosities ranged from 0.200 to 1.000 and 0.683 to 0.879, respectively. We tested for cross-amplification of these 10 markers in the related species L. hirsuta and found that nine could be successfully amplified. CONCLUSIONS: The SSR markers characterized here are the first to be developed and tested in L. tibetica. They will be useful for future population genetic studies on L. tibetica and closely related species.